Draft version February 5, 2008 

Preprint typeset using I4TgX style emulateapj v. 9/08/03 



THE FLAMINGOS EXTRAGALACTIC SURVEY 

Richard J. Elston'-^, Anthony H. Gonzalez'-\ Eric McKenzie', Mark Brodwin*, Michael J. I. Brown' '', 
Gustavo Cardona', Arjun Dey\ Mark Dickinson', Peter R. Eisenhardt'*, Buell T. Jannuzi', Yen-Ting Lin', Joseph 
J. Mohr'-*, S. Nicholas Raines', S. A. Stanford'"', Daniel Stern* 

Draft version February 5, 2008 



in 
o 
o 

(N 

> 
O 



> 

(N 



in 
o 



Oh 
I 

O 



X 



ABSTRACT 

Using the Florida Multi-object Imaging Near-IR grism Observational Spectrometer (FLAMINGOS), we have 
conducted the FLAMINGOS Extragalactic Survey (FLAMEX), a deep imaging survey covering 7.1 deg^ 
within the 18.6 deg^ NOAO Deep Wide-Field Survey (NDWFS) regions. FLAMEX is the first deep, wide- 
area near-infrared survey to image in both the J and K., filters, and is larger than any previous NIR surveys of 
comparable depth. The intent of FLAMEX is to facilitate the study of galaxy and galaxy cluster evolution at 
1 < z < 2 by providing rest-frame optical photometry for the massive galaxy population at this epoch. This 
effort is designed to yield a public data set that complements and augments the suite of existing surveys in the 
NDWFS fields. We present an overview of FLAMEX and initial results based upon ^ 150,000 TSQ-selected 
sources in the Bootes field. We describe the observations and reductions, quantify the data quality, and verify 
that the number counts are consistent with results from previous surveys. Finally, we comment upon the utility 
of this sample for detailed study of the ERO population, and present one of the first spectroscopically confirmed 
z > 1 galaxy clusters detected using the joint FLAMEX, NDWFS, and Spitzer IRAC Shallow Survey data sets. 

Subject headings: surveys — infrared: galaxies — galaxies: evolution — galaxies: clusters 



L INTRODUCTION 

Two fundamental goals of extragalactic astronomy are to 
understand the star formation and mass assembly histories of 
the universe. It has become increasingly apparent in the last 
few years that the epoch z= 1 - 2 is critical for the evolution 
of both quantities. Current observations indicate that the bulk 
of stellar mass assembly in galaxies and ^ 50% of the to- 
tal star formation occurs in thi s interval (Fo ntana et al. 2003j 
[Dickinson et al. 2003; Rudni ck et alJl2003l) . It is also dur- 
ing this epoch that we expect to witness the formation of the 
first massive galaxy clusters (for example, see|Haiman et al. 
llpOl). Observations at this epoch thus promise to provide a 
portrait of the universe during its most vibrant evolutionary 
period. Yet, in c onstrast t o the wealth of data at low redshifts 
\z<l, 'York et al. 2002^ iCoUess et al."2001f: 'Kochan ek etall 
1^04) and high redshifts (z = 2.5-4, Steidel et al. 1998. 119991: 
[Adelberger et al. 2003*), this redshift regime remains less well 
explored due largely to the observational challenges histor- 
ically associated with the "redshift deser t" at z = 1.4-2.5 
(for an overview, see 'Steidel et al."2004^. Recent s pectro- 
scopi c work by a number of groups (e.g., Abraha ni et alJ 
12004; Coil et al. '2004; Steidel et al. 2004; Adelberger et aj 
12005 a b; Mignoli et al. 2005) has begun to open this redshift 
window, but we currently lack large, homogenous samples 
comparable to those that exist at other redshifts. Moreover, 
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these existing programs are restricted to relatively small ar- 
eas and are thus unable to probe the fraction of the massive 
galaxy population that resides in the overdense group and 
cluster environments. Infrared imaging surveys are a key first 
step towards generation of such samples. Infrared observa- 
tions (both near-infrared and with IRACon the Spitzer Space 
Telescope) have the advantage that the combined evolutionary 
and A:-corrections are mild at this epoch, which means that 
uniform, roughly stellar mass-limited samples can be selected 
over a range of redshifts. 

Using the Florida Multi-object Imaging Near-IR grism 
Observational Spectrometer (FLAMINGOS), we have con- 
ducted the FLAMINGOS Extragalactic Survey (FLAMEX). 
FLAMEX is a wide-area, near-infrared (NIR) imaging sur- 
vey that is motivated by two basic considerations. First, 
FLAMEX is designed to probe the evolution and cluster- 
ing of the massive galaxy population at z = 1-2 and to 
identify a sample of galaxy clusters at this epoch that can 
be used to study the assembly of cluster galaxies. The 
NDWFS fields are the targets of one of the most ex- 
tensive panchromatic investigations in all astronomy, with 
space- and ground-based imaging and spectro scopic programs 
spanning radio to X-ray wavelengt hs (Jannuzi &DeiJn992 
iRhoads et al."2000"; de Vries et al. 2002; Hoopes et al.'"200a 
Lonsdale et al. 2003; Eisenhardt et al. 2004; Kochanek et all 
2004; Pierre et al. 2004; iHouck et alJ 1200^: IMurrav etall 
i2005) . The northern (Bootes) field in fact is the only wide- 
area survey region presently mapped by Chandra, GALEX, 
and Spitzer (with IRAC and MIPS). 

FLAMEX provides a combination of area and depth in both 
J and Ks that is unique among current NIR surveys (see Table 
Hand Figure Coupling the FLAMEX data with the opti- 
cal (and when available IRAC) photometry enables derivation 
of robust photometric redshifts (Brodwin et al., in prep., and 
Brown et al., in prep.), providing an efficient means of iso- 
lating large samples of massive galaxies at z = 1 - 2. This 
combination of mass selection and photometric redshifts over 
a wide area constitutes a powerful tool not only for galaxy 
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Table 1. NIR Surveys 



Survey 



Skrutskie et al. f 1997, 2MASS) 

Gardner et al. (1996) 

Jannuzi & Dev (199^ NDWFS Bootes) 

Huana et al. (1997^ 

FLAMEX 

Chen et al. (2002. LCIRS) 
Drorv et al. (2001, MUNICS) 
Kummel & Waaner (200r) 

Szokolv et al. (1998) 
Huang etal. (2001. CADIS) 
Daddi et al. (2000) 

Glazebrook et al. (1994) 
Cristobal-Homillos etal. (2003) 
Minezaki etal. (1998) 

Martini (2001) 
^imatti et al. (2002, K20) 
McLeod et al. (1995^ 

Totani et al. (20011 
piorgovski et al. (1995) 
Moustakas etal. (1997) 
^ershadv et al. (1998) 



Filter 


Area (deg^) 


Depth 


JHK, 


All Sky 


K=U.5 


K 


9.84 


K= 15.6 


K 


9.3 


K= 18.6 


K' 


8.23 


K' = 16 


JKs 


7.1 


Ks f» 19.3 


H 


1.1 


U ^ T\ C 


J IS. 


~1 


' — 1 Q S 
iv — ly.J 


K 


0.91 


K=n 




0.61 


K= 17.5 




0.6 


Ks = 16.5 


K' 


0.2 


K' = 19.75 




0.19 


^Cs = 18.8 




0.12 


K= 19.2 


K 


0.15 


17.3 


K, 


5 X 10-^ 


Ks = 2l 


K' 


5 X 10-2 


K' = 19 




6x 10^ 


= 21 


JHK 


3.4 X 10"- 


/(■= 18.5 




1.4 X 10"- 


/f >20 


K 


5.6 X 10-3 


K = 20 




5.6 X 10-^ 


/f = 21.5 


K' 


1.1 X \0-^ 


^r' = 23.5 


K 


8.3 X 10^ 


^: = 24 


K 


5.6 X 10^ 


K Si 23.4" 


JK 


4.2 X 10-* 


K = 2A 



Note. — This table is intended to be a representative rather than comprehensive 
synopsis of previous surveys. The surveys are listed in order of decreasing survey area. 
We note that different authors quote different signiticance limiting magnitudes (for ex- 
ample 3(7 vs. 5(7) within different size apertures, so the listed values provide only an 
approximate comparison of survey depths. 

"This survey is split between two regions with limiting magnitudes of K = 
23.1 and A' = 23.7. 



cluster detection and galaxy evolution studies, but also for se- 
lection of subsamples of rare objects that can be studied in 
greater detail. 

Second, the addition of FLAMEX to the NDWFS reper- 
toire provides the requisite deep NIR photometry that will 
be required by the next generation of NIR spectrographs. 
The astronomical community is poised to enter a new era in 
NIR spectroscopy, with multiobject instruments coming on 
line at multiple 8-lOm class telescopes (MOIRCS at Subaru, 
FLAMINGOS II at Gemini South, EMIR at the Gran Telesco- 
pio Canarias). In this regard, the FLAMEX survey furthers 
one of the original aims of the NOAO Deep Wide-Field Sur- 
vey by providing a wide-area public data set that can be used 
to select targets for spectroscopy. 

In this paper we present an overview of the FLAMINGOS 
Extragalactic Survey. The science programs based upon the 
survey and details of the FLAMEX catalogs will be presented 
in other papers. The main aim in this manuscript is to provide 
relevant information regarding the general survey character- 
istics to the community in conjunction with the public data 
release, enabling researchers interested in the NDWFS fields 
to assess the applicability of this data set to their research pro- 
grams. Towards this end, after presenting the data reductions 
and data quality we use the catalog in the Bootes field to pro- 
vide two brief illustrations of the science topics that can be 
addressed with the FLAMEX data set. 

The layout of this paper is as follows. We begin by describ- 
ing the experimental design and observations in ^ including 
target sensitivities and field geometry. In ^we present the 
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Fig. 1 . — Comparison of FLAMEX with previous K-hand surveys. The 
points show the published limiting magnitudes and areas for all the surveys 
listed in Table [TI except for the //-band LCIRS. Among the wide-area sur- 
veys shown, only FLAMEX, 2MASS, and MUNICS provide complete cov- 
erage in both J and Ks data. We note that no coirection has been applied to 
correct for the different filters (K, K' ,Ks) used by the different surveys, and 
that different authors quote different significance limiting magnitudes, so the 
data points provide only an approximate comparison of survey depths. 



reduction procedure, followed by a description of the proper- 
ties of the final Bootes catalog in ^ Illustrations of initial 
science and a discussion of the data release are given in ^ 
and the paper is summarized in ^ When relevant we assume 
a standard cosmology with il^f = 0.27, J7a = 0.73, and Hq = 70 
km s"'Mpc"'. 

2. EXPERIMENTAL DESIGN AND OBSERVATIONS 

This survey is made possible by the FLorida Multi-object 
Imaging Near-Infrared Grism Observational Spectrometer 
(FLAMINGOS), which when commissioned was the world's 
first cryogenically cooled multi-object infrared spectrograph 
('Elstonet al."2003y FLAMINGOS has a fast all-refractive 
optical system that can be used at telescopes slower than f/7.5. 
This optical system, coupled with a 2048x2048 HgCdTe 
HAWAII-2 array, makes FLAMINGOS a very efficient wide- 
field imager when used on fast, small aperture telescopes. For 
the FLAMEX survey we rely purely on FLAMINGOS' imag- 
ing capabilities and the wide field afforded by the Kitt Peak 
National Observatory 2.1m. At this telescope FLAMINGOS 
has a 21' X 21' field-of-view with 0.61" pixels. 

2.1. Location and Field Geometry 

We initially targeted rectangular 5 deg^ regions within each 
of the two NDWFS fields. The FLAMEX survey chose the 
NDWFS fields due to the deep optical imaging that was in 
progress in these fields when the survey was begun, with the 
BwRI data being essential for enabling derivation of accurate 
photometric redshifts. The even split between the two fields 
was designed to permit follow-up observations during both 
observing semesters and from both hemispheres, while still 
yielding coverage for large, contiguous regions. Finally, the 
angular size of the survey is driven by the requirement that we 
probe a sufficiently large volume to detect massive clusters 
at z > 1. For a standard ACDM cosmology, one expects a 
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Fig. 3.— Final FLAMEX coverage in the NDWFS Cetus field. The 
image and boxes are the same as in Figure l2l except that three fields with 
incomplete /-band coverage are drawn in white. For Cetus we also overlay 
the locations of the SWIRE survey (large, white box), LCIRS (large, thin 
black box), and GDDS (small, thin black box). Near the GDDS field the 
overlap of the FLAMEX pointings was increased to increase the depth. The 
survey covers a total of 2.4 deg^ in both J and K, in Cetus, with an additional 
0.3 deg2 in K,. 



Fig. 2. — Final coverage for the FLAMEX survey in the ND- 
WFS Bootes field, overlaid upon an IRAS 100/im image of the region 
pVheelock et al. 1994). The boxes denote individual FLAMEX pointings, 
with the extent of the full NDWFS field also displayed. In the Bootes region 
we cover a 4.7 deg^ strip. 



surface density of approximately 0.7-4 deg"^ clusters with 
M > 10'^ Mq at z > 1, depending upon the value of erg." 
Consequently, the survey must cover roughly 5-10 deg^ to 
ensure detection of a sample of massive clusters at this epoch. 

While the survey plan called for a total of 10 deg^, in prac- 
tice worse than expected weather (see below) limited the total 
coverage to 7.1 deg^. In the northern (Bootes) field the sur- 
vey covers a 4.7 deg^ region in both J and K,; in the southern 
field (Cetus) the corresponding coverage is 2.4 deg^. The fi- 
nal geometry of each field is shown in Figures |2j]2[ In Cetus, 
this geometry was driven by following considerations. First, 
we started with observations in the eastern part of the ND- 
WFS field to best complement the shallower j T-ban d imag- 
ing from the NDWFS survey (Jannuzi & Devl lT999l Dey et 
al., in prep.), which predominantly lies in the western part 
of the Cetus NDWFS field. Second, we attempted to both 
maximize overlap with the Spitzer Wide-Area Infrared Extra- 
galactic Survey (JLonsdal e et al. 2003, SWIRE), which over- 
laps with the eastern part of the NDWFS, and cover the Gem- 
ini Deep Deep Survey's ("Abra ham et al..2 004. GDDS) field in 
Cetus. The former provides valuable additional information 
for photometric redshifts and constraining galaxy populations 
at z > 1, while the latter provides calibration for photometric 
redshifts in the redshift desert (z w 1 .5-2). Observations were 
conducted as a gridded series of pointings, with adjacent field 
centers separated by 20'. 

2.2. Survey Specifications 

The original technical specifications for the program pre- 
scribed imaging of the survey region in both J and Ks for 
2 hours to reach depths of 7 = 22 and K, = 20.5 (5ct, Vega) 
within 2" diameter apertures, assuming 1.5" FWHM seeing, 

" These surface densities are calculated a nalytically, and are insen si- 
tive to whether one uses the prescription from ISheth & TormenI 11999ft or 
IJenkins etalHSOOih . 



which for a point source is equivalent to total magnitudes 
7 = 2 1.4 and Ks = 19.9. The final catalog is shallower (see 
^4.1> due to both the instrumental aberrations and worse than 
expected image quality. The effective depth of the data are of 
course also a strong function of the ambient temperature (i.e. 
the sky brightness). We therefore initially imaged each field 
for two hours per filter, and then reimaged the shallower fields 
in later runs to obtain relatively uniform data. 

2.3. Observing Procedure 

The standard observing sequence for the survey consisted 
of a randomly ordered 5x5 dither pattern, with 1' offsets 
about the pointing center. This pattern was repeated until the 
required total exposure time was achieved. Individual expo- 
sures times were dependent upon the ambient temperature, 
ranging from 20-60s in and 90- 120s in J. The typical 
background sky level in Kg was 20,000-25,000 counts, with a 
maximum allowed level of ^^30,000 counts. Detector nonlin- 
earity can be corrected to better than 1% up to 45,000 counts; 
the maximum allowed level is designed to ensure that 2MASS 
stars with Ks > 12.5, which are used for photometric calibra- 
tion, have peak fluxes below this level. 

2.4. Observing Conditions 

Data were obtained on the KPNO 2.1m during an alloca- 
tion of 97 nights spanning 2001 December to 2004 Decem- 
ber While photometric conditions were not required, us- 
able data were obtained on only ^ 50% of these nights due 
to a combination of suboptimal weather conditions and tele- 
scope mechanical failures. Observations were coordinated 
with the NOAO Survey Program "Toward a Complete Near- 
Infrared Spectroscopic and Imaging Survey of Giant Molec- 
ular Clouds" (PI: Lada). The two surveys shared nights to 
maximize observing efficiency since the target fields for the 
two groups lie at complementary right ascensions. 

3. REDUCTIONS 

Survey imaging was reduced using the LONGLEGS 
pipeline (Roman et al., in prep., and Gonzalez et al., in prep.), 
which employs standard infrared reduction procedures. Lin- 
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Fig. 4. — Map of the 50% completeness limit for 6" aperture magnitudes 
for Ks in one of the Bootes subfields. North is down and east is to the left. 
This map is generated using simulated elliptical galaxies with sizes consisten t 
with the observed galaxy size distribution at z = 1 - 2 I Ferguson et all2004l) . 
As discussed in the text, the survey depth can vary by ~ 0.5 magnitudes 
within a given subfield. 



earity correction, dark subtraction, and flatfielding were all 
performed in the normal fashion. A running sky, typically 
comprised of 8 adjacent frames, was used for sky subtraction, 
with objects detected in the individual images masked before 
creation of the sky frames. The images were then transformed 
to correct for geometric distortion, doubling the number of 
pixels along both axes, and combined using a integer shift- 
and-add approach on these oversampled images. The result- 
ing stacks are astrometrically and photometrically calibrated 
using th e Two Micro n All Sky Survey (2MASS) All Sky Data 
Release (iCutri et al.ll2b03D . The final photometric calibration 
was computed using the average photometric offset for color- 
selected 2MASS stars in each field, weighted by the photo- 
metric errors. Field-to-field variations in the J-K^ color of 
the full stellar locus are typically at the level of 0.02 magni- 
tudes or less, indicating that the calibrations for the individual 
filters are good to this level or better 

Because the data for the survey were obtained over an ex- 
tended time period, in many instances it was necessary to 
combine image stacks from different nights (or years). When 
combining, we weighted the images based upon the seeing, 
zeropoint, and rms noise in the input frames via the formula: 

]^Q0.4mo 

aiFWHMf 

where mo is the photometric zeropoint, FWHM is the full 
width half-maximum of the point spread function (PSF) at the 
field center, and a is the rms noise in each image. 

Details related to the FLAMINGOS linearity correction and geome tric 
distortions can be found at |http : / / flamingos ■ astro . ufl ■ edu| 




Completeness Limit (6" aperture) 



Fig. 5. — Survey area in Bootes as a function of the 50% completeness 
limit for the 6" diameter aperture. Over 90% of the survey region is complete 
to Ks = 19.2, with 50% being complete to Ks = 19.5. Note that while the total 
area covered in Bootes is 4.7 deg^, the effective ai'ea after accounting for 
masking and holes in the survey coverage is 4.55 deg^ . 

4. CATALOG PROPERTIES 

We construct /Tj-selected catalogs for each survey sub- 
field using So urce E xtractor version 2.3 (SExtractor, 
Hertin & Arnout3 fT996l) . using dual image mode to measure 
the y-band photometry within the same regions. Detection 
is performed in Ks with a 0.76" FWHM gaussian convolu- 
tion kernel, and we impose a 5a object detection threshold. 
Weight maps are also used to minimize detection of spurious. 
The catalogs for the individual subfields are then merged, with 
strict right ascension and declination boundaries defined be- 
tween subfields to avoid double-counting of objects. We use 
completeness simulations, described below, to compute more 
realistic estimates of the magnitude errors than those provided 
by SExtractor, since the simulations can coiTectly account for 
correlated noise across pixels and the impact of background 
fluctuations. The typical magnitude errors from the simula- 
tions are roughly 2-3 times as large as the SExtractor errors 
for most apertures, and larger for the 20" aperture and the 
SExtactor automatic aperture ("AUTO"), which is designed 
to give an estimate of the total magnitude. 

4. 1 . Survey Depth 

The depth of the survey is position-dependent due to both 
instrumental aberration and shorter effective exposure times 
near the stack edges. The most appropriate way to quantify 
this variation is via completeness simulations for mock cat- 
alogs. We perform such simulations for each field to gener- 
ate 50% completeness maps. Schematically, the design of the 
simulations is as follows. 

First, we require a two-dimensional model for the point 
spread function for each field. For this we take an empirical 
approach - we use data from the dense stellar fields observed 
by the star formation survey (see ^2.4t during the same ob- 
serving run. We select the individual frames with the best 
seeing, and from these data generate a lookup table of the co- 
ordinates of bright, unsaturated stars. Several of the runs lack 
imaging from the star formation survey; for these we boot- 
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Fig. 6. — Infrared color-magnitude diagram for tlie Bootes field. The 
figure includes all objects with Ks < 20.5 that do not overlap with other 
objects (~84,000). Objects are grouped in bins of width AKs = 0.1 mag 
and AJ -Ks = 0.035 mag. The greyscale runs from 0-100 objects per bin, 
with a logarithmic intensity scale. Note that the two horizontal sequences 
at J — Ks = 0.4 and J — Ks = 0.8 con'espond to Galactic sta rs with types later 
than G5 and earlier than K5, respectively i Finlator et all200Cl) . Intermediate 
type stars (K5 through G5) lie between these two sequen ces. The two curves 
are th e predicted magnitudes and colors, based upon Bruzual & CharloJ 
120031) models with the Padova 1994 isochrones ( Bertelli et al. 1994), for 
M = 2 X lO" M0 galaxies with single burst (solid) and exponentially de- 
clining (dashed, r = 10 Gyr) star formation histories. The points correspond 
toz = 0.25,0.5,l,1.5,2. 



Strap off the survey data, using stars from the survey fields 
with the best seeing. To test the reHabihty of this bootstrap 
approach, we also run bootstrapped simulations for several 
fields with stellar template data. In these instances the com- 
pleteness limits agree to < 0.05 mag in the field center and 
< 0.1 -0.15 mag in the corners. 

Second, we randomly select a location in the final image to 
insert a fake galaxy, and compute the corresponding location 
in each individual frame. We then generate a fake galaxy, con- 
volve it with the nearest PSF star for each individual frame, 
combine these convolved model images, and add the galaxy 
into the final image. 

For every field, in both the J and filters, we insert 
5000 fake galaxies with magnitudes spanning the complete- 
ness limit. The effective radii of these galaxies are chosen 
to uniformly span th e range of values seen at z = 1 - 2 by 
IFerguson et alJ J2004) . who measured the distribution of half- 
light radii of distant galaxies in the Hubble Deep Field-North 
and Chandra Deep Field-South. Because the primary science 
goals of the survey focus upon the elliptical galaxy population 
at z > 1 , we do not attempt to reproduce the full morphologi- 
cal distribution seen at this epoch. Instead, all mock galaxies 
used to calculate the completeness limits are pure bulges, with 
a distribution of effective radii in the range 0.07-0.5" (0.5- 
4 kpc at z = 1). While obviously a simplification, in practice 
the bulk of the objects detected in the survey are unresolved in 
the FLAMEX imaging data, making this a reasonable approx- 
imation even for disk galaxies. To verify this assertion, we 
reran the simulations in two subfields using only spiral galax- 
ies with scale lengths = 0.6" (4.8 kpc at z = 1). The derived 
completeness limits were consistent to within 0.1 magnitudes 
with the results for pure bulges. 
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Fig. 7. — Color-color diagram for all FLAMEX sources, using 6" diameter 
aperture magnitudes. The solid line denotes the division between stars and 
galaxies that we utilize in il4.3l 



We use the output of the simulations to derive two- 
dimensional maps of the 50% completeness limit for aperture 
and AUTO magnitudes measured with SExtractor. Figure |4] 
shows the derived map for a representative subfield. As ex- 
pected, the off-axis aberrations degrade the survey depth in 
the corner regions. Statistical uncertainties for these maps 
can be derived via bootstrap resampling; in the image centers 
the completeness values have uncertainties of < 0.03 mag, 
while in the corners the uncertainties are '^0.1 mag. These 
values are slightly smaller that the systematic uncertainties 
evidenced by comparison of simulations using the bootstrap 
stellar templates with those from the star formation fields (see 
above). 

In Figure [S] we plot the survey area as a function of the 
50% completeness depth in the Bootes region for 6" diameter 
aperture magnitudes, incorporating the exposure maps to re- 
produce the effect of masking and dithering. More than 90% 
of the survey region is complete to Ks = 19.2, with 50% being 
complete Xo Ks = 19.5. 

4.2. Star-Galaxy Separation in Color Space 

FLAMEX is the first deep, wide-area near-infrared survey 
conducted in two filters. Our primary motivation for obtaining 
both / and data was to efficiently isolate the z > 1 popu- 
lation, and these two bands, when combined with the optical 
NDWFS data (and/or with IRAC data), are particularly valu- 
able for deriving robust photometric redshifts at z > 1 (Brod- 
win et al., in prep.; Brown et al., in prep.). The combination of 
near-infrared (NIR) and optical photometry from FLAMEX 
and the NDWFS also find a number of other uses, including 
star-galaxy separation. 

The variable PSF noted above precluded use of structural 
information to separate the star and galaxy populations in the 
FLAMEX data, but color information is an effective substi- 
tute. The J-Ks color alone provides a simple means of per- 
forming star-galaxy separation, as can be seen in the color- 
magnitude diagram (Figure|6|l. Galactic stars with types ear- 
lier than K5 lie at J-Ks = 0.8; later types have bluer colors, 
with types later tha n M5 occupying a locus al J-Ks = 0.4 
dFinlator et alJl2000h . Galaxies, which occupy the redder lo- 
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Fig. 8. — Ks differential number counts for the Bootes region. We use 
Source Extractor AUTO magnitudes, wliicli for blight extended sources agree 
well with the 2MASS K20 total magnitudes. The thin solid line denotes the 
total number counts in the region, while the thick solid and dotted lines cor- 
respond to the color selected galaxies and stai's, respectively. We also include 
the completeness-corrected galaxy number counts (dot-dashed line) and the 
galaxy n umber coun ts from the 20 arcmin^ Calar Alto Deep Imaging Survey 
(CADIS ;'Huang et al. 2001), one of the largest previous Kg surveys to compa- 
rable depth. Our completeness-corrected number counts are consistent with 
the CADIS values over the full range, albeit with greater statistical precision 
(see also Figurel9l. 



cus of objects, are well-separated from the stars for z > 0.05. 

While the J-K^ color alone provides a reasonable means of 
performing star-galaxy separation, inclusion of the NDWFS 
optical photometry can be used to improve the separation. 
We have generated a matched catalog combining the NDWFS 
third data release (DR3; Jannuzi et. al., in prep) and FLAMEX 
data. Detections in the different bands were matched if the 
centroids were within 1" of each other For extended objects, 
detections in the different bands were matched if the centroid 
in one band was within an ellipse defined using the second 
order moments of the light distribution of the object in an- 
other band.'^ FigureQshows the Bw-I vs J-K, color-color 
diagram for the survey using 6" aperture magnitudes (suffi- 
ciently large that PSF mismatch between surveys is not an 
issue). The stellar sequence can be seen to cleanly separate 
from the galaxy locus, with the M stars being relatively redder 
in Bw - 1 than the bluest galaxies, and the G stars remaining 
bluer in J- Ks. There remains a modest overlap between faint, 
nearby galaxies and the G star population, but the cut shown 
in Figure overall provides an efficient means of separation 
w ith < 5% misidentification of galaxies at faint magnitudes 

4.3. Number Counts 

Galaxy number counts provide a useful means of te sting 
the reliability of both the star-galaxy separation ( ^4.21 and 
the completeness simulations ( 34. H . In Figure|S]we show the 
differential number counts for all sources in Bootes. Using 
the star-galaxy color separation, we split the number counts 
into the star and galaxy contributions. The slight upturn in the 

This ellipse was defined with the SExtractor parameters 2 X 
A_WORLD, 2 X B_WORLD, and THETA_WORLD. 



Stellar counts at A!^, w 19 indicates that we may be misclas- 
sifying up to 5% of the galaxy population as stars if this is 
not an intrinsic feature of the stellar number counts. This sys- 
tematic should be considered for any work involving either 
high-precision modelling of the number counts or statistical 
analysis of the local faint galaxy population. 

Figure|8]also provides a comparison of the FLAMEX num- 
ber counts with results from the 0.2 deg^Calar Alto Deep 
Imaging Survey (CADIS, Huang et al. 2001), one of the 
largest previous Ks surveys to comparable depth. Once a 
completeness correction is applied to the galaxy counts (dot- 
dashed curve), the faint-end FLAMEX number counts are 
consistent with CADIS to Ks > 19.5. This agreement at the 
faint end provides validation for the completeness simulations 
described in 34.11 

Near-infrared galaxy number counts are also a classic di- 
agnost ic for constraining models of galaxy evolution (for ex- 
ample, Moustakas et al. ■n^. "Hua ng et alJl200U IChen et alJ 
|2£)02; Cristobal-Hornillos et al. 2002l)- At these wavelengths 
the K-corrections are small, and the luminosities are both less 
affected by dust and less sensitive to recent star formation 
than optical number counts. These factors result in samples 
that are closer to mass-selected and hence simpler to model. 
While the magnitude regime probed by the FLAMEX data 
is not new, the large area afforded by this survey provides 
significantly greater statistical precision than previous work, 
enabling refined modelling. Realistic semianalytic models, 
which attempt to simultaneously match the observed counts 
in Ks and other passbands (such as Nagashima et al. 2002), 
are the optimal means for using this information. Such a de- 
tailed analysis is beyond the scope of this paper; however, in 
Figure |9] we compare with a simple classical model to illus- 
trate both the utility of this data set and its precision relative 
to previous studies. 

We use extendable galaxy number count model, ncmod, 
from Gardner (1998) with the same spectral templ ates and 
morphological mix of ellipti cals and spirals as iGardned 
(il99&) . Similar to .Cristobal-Hornillos et alJ J2003h . we take 
the formation redshifts to be z = 2 for ellipticals and z = 1 
for spirals, and add a component of star-forming dwarf galax- 
ies with faint end slope a = -1.5. For the general galaxy 
popul ation we u se the Schechter function parameters from 
Koch anek et al.l(l2001.) . while for the star-forming dwarf pop- 
ulatio n we use the same prescri ption as Gardner (1998) and 
iCristobal-Hornillos et al.](l2003h . but with = 3.5 x lO'^h]^ 
Mpc"^ (roughly a factor of two lower than these studies). As 
found by Cristobal-Hornillos et al. (2003), this type of model 
qualitatively fits the number counts over much of the observed 
magnitude range, albeit with some residual deviation. Most 
notably, FLAMEX and the other surveys yield fewer counts 
than the model at the brighest magnitudes. The FLAMEX 
number counts also exhibit an excess at the faint end {Ks > 19) 
that is statistically significant with regard to the plotted Pois- 
son uncertainties, albeit at a level that is small compared to the 
statistical uncertainties for most of the other samples shown. 
While an excess is also evident in the faintest CADIS data 
point, we caution against overinterpretation of this feature be- 
cause this upturn occurs in the regime where our completeness 
correction is significant. 

5. DISCUSSION 

The fundamental motivation for the FLAMINGOS extra- 
galactic survey is to facilitate systematic study of the evo- 
lution of massive galaxies and galaxy clusters at z = 1-2, 
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Fig. 9. — Differential number counts from FLAMEX and a representative sample of published data sets fGardner et alj |1993t IDiorgovski et 311119951 
IGardner et al. 1996; Huane et al. 1997; Moustakas et al. 1997; Bershadv et al. 1998; Huang et al. 2001 ; Kochanek et al. 2001 ). The left panel shows the number 
counts for the full magnitude range probed by existing studies. The FLAMEX data and a subset of the data sets are shown in the upper part of the figure. Overlaid 
is a simple model constructed using ncmod that is similar to the one tested in Cristobal-Hornillos et al. (2003). The formation epoch for ellipticals and spirals 
are set at z = 2 and z = 1, respectively, and a population of star-forming dwarfs (with a lower 0, = 2.8 X Kr'/i^ Mpc"') is included to reproduce the faint end 
number counts seen in deep surveys. The lower part of the left panel shows the residuals relative to this model, including the en'or bai's, for all data sets. For 
FLAMEX the error bars include only the Poisson uncertainty. There is good agreement between FLAMEX and previous studies in the magnitude range probed 
by FLAMEX. To facihtate a clearer comparison, the panels on the right show the same differential number counts, but restricted to the magnitude regime probed 
by FLAMEX. The different surveys are generally consistent to within the statistical uncertainties. 



the epoch of peak global star formation and a critical era 
for both galaxy and galaxy cluster assembly. Complement- 
ing the existing programs in the NDWFS fields at other wave- 
lengths, FLAMEX provides rest-frame optical photometry for 
the massive galaxy population at this epoch. The combina- 
tion of the FLAMEX and NDWFS data enables the robust 
photometric redshifts and uniform selection criteria necessary 
for study of the galaxy and galaxy cluster populations at this 
epoch. In Bootes, the FLAMEX data also serves to bridge the 
gap between the existing optical and IRAC data sets, facili- 
tating cross-association and investigation of optical drop-out 
populations. 

The key scientific issues listed above will be the focus of a 
series of separate papers. In the following paragraphs, how- 
ever, we supplement this overview of the survey with a few 
illustrations of basic applications of the FLAMEX data to pro- 
vide the reader with a more quantitative sense of the science 
that can be achieved. In particular, we provide a brief analysis 
of the extremely red object (ERO) population and a synop- 
sis of our galaxy cluster search before discussing our public 
data release. These items are intended simply to highlight a 
few uses of the FLAMEX data set - examples of other, mul- 
tiwavelength programs (either underway or planned) that use 
FLAMEX include a high-redshift quasar search, investigation 
of the morphology-density relation at z > 1, construction of 
empirical optical to mid-IR galaxy spectral energy distribu- 
tions, and investigation of the ensemble properties of distant 
galaxy clusters. 

5.1. EROs 

First discovered bv lElston et alJ il988l) . extremely red ob- 
jects (EROs) are traditionally defined as objects W\i\\R-Ks > 
5 (althoug h a range of definitions exist in the literature; see 
lMcCarthvll2()Ql . While initially speculated to correspond 



to a population of protogalaxies at z > 6, the ERO popu- 
lation is now understood to be predominantly comprised of 
two distinct types of galaxies - dusty starbursts and galax- 
ies dom inated by old stellar populations, both at z ~ 0.8- 2 
(e.g., Elston et al. 1989; Cimatti et al. 2002; Yan et al. 2004«). 
For the dusty starbursts, strong extinction in the rest-frame 
near-ultraviolet drives the observed color; for the old galax- 
ies, the red color is a simple consequence of the 4000 A break 
passing redward of /?-band. These two populations together 
are believed to constitute the dominant progenitor population 
for present-day elliptical galaxies, a belief that is supported 
by their space densities jP addi et al. 2000) and strong spa- 
tial clustering ('paddiitdJ|200Qf |Brown.eL^ but see 
Yan et al. 2004 for concerns with this interpretation). How- 
ever, as a consequence of the strong spatial clustering and rel- 
atively low surface density of EROs, large area surveys are 
required for precise statistical analysis of the population. 

The FLAMEX survey provides an order of magnitude in- 
crease in EROs sample size compared to other recent studies. 
We define EROs as being objects with R-Ks > 5 within a 6" 
aperture, imposing the additional requirement that J-Ks> 1.2 
to eliminate potential contamination from stellar sources with 
spurious optical photometry (Figure llOb). By this defini- 
tion, in the Bootes field alone we have over 7,000 EROs at 
Ks < 19.5 (Figure [Tot'); for comparison, the largest sample 
used thus far to study the clustering of EROs contained 67 1 
objects down to K =18.4 within a 1.4 deg^ region in the ND- 
WFS Bootes field (iBrown et alJ2005 ). The FLAMEX counts 
are consistent with previous work, with the advantage that the 
large field enables identification of the brightest, rarest EROs 
and measurement of the number counts as bright at A'j w 17. 
This sample is also sufficiently large to begin to quantify the 
number density distribution as a joint function of magnitude 
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Fig. 10. — ((a) R-Ks versus J—Ks color distribution of FLAMEX sources in Bootes for a subset of sources with 17 < ^f, < 19. EROs are identified as objects 
with J — Ks > 1.2 and R — Ks > 5, with all colors calculated using 6" aperture magnitudes. Filled triangles in this figure con'espond to objects with R > 24.5 
(roughly equivalent to an i?— band SNR=4 in a 6" aperture), and should be considered lower limits on the R - K, color. Roughly 40% of the EROS have lower 
limits m R-Ks. Open boxes correspond to objects that qualify as EROs from their R — Ks color, but have only lower limits on the J -Kg color The vertical 
dashed line indicates that 50% of sources with J—Ks > 2.4 are undetected in the /—band magnitude, (fo) Cumulative surface density of EROs detected within a 
4.2 deg^ subregion of the Bootes field. For R—Ks > 5 we present the FLAMEX Bo otes data as filled circles and compare the observed number density with the 
results of several previous sfitdies, including'Paddi et al.'>2O0lll),IR oche er^l200l , and'B rown et"ai]l200 l. The FLAMEX counts for redder cuts (R-K, > 5.5 
and R — Ks > 6) are presented as solid lines down to limiting magnitudes corresponding to R = 24.5. We note that these curves are both roughly parallel to the 
R—Ks>5 data. The total number of EROs in the Bootes field is shown on the right axis, demonstrating that the sample includes > 7,000 EROs with Ks < 19.5. 
We apply no completeness coirection to the FLAMEX data, and hence the true suii'ace density exceeds the FLAMEX points for A^, > 19 (see Figure|8). There is 
also no correction applied to convert some of the other surveys from KtoKs, but for EROs the difference is minimal. 



and color (see the curves for redder cuts in Figure [TOJ?), pro- 
viding a more stringent observational test for semianalytic 
evolutionary models. 

5.2. Galaxy Clusters 

One of the principal motivations for FLAMEX is identifi- 
cation of a large sample of galaxy clusters at z > 1. In the 
standard paradigm of hierarchical structure formation, we ex- 
pect the formation of the first massive galaxy clusters to oc- 
cur at at z = 1-2. The comoving cluster number density at 
M > lO'"^ M0 increases by roughly a factor of 30 between 
z = 2 and z = 1 (WMAP cosmology), with the precise num- 
ber density and degree of evolution strongly sensitive to the 
details of the underlying cosmological model (particularly ag 
and r^/u) - Aside from the cosmological applicati ons (see for 

example iHaiman et al.ll200ll iBor^ani eT al 200 li iWan^ et al.l 

l2004h . detailed studies of homogeneous cluster samples at 
this epoch are essential for understanding the evolution of the 
cluster galaxy population and how the interplay between the 
galaxies and host clusters impacts this evolution. Observa- 
tions at this epoch constrain the evolution of the halo occu- 
pation distribution (HOD), the formation epoch of cluster el- 
lipticals, evolution of the color-magnitude relation, and the 
morphological evolution of cluster galaxies. 

Infrared imaging, coupled with shorter wavelength photom- 
etry, provides an efficient means of detecting clusters at this 
epoch. The infrared data enable construction of roughly stel- 
lar mass-selected galaxy samples, and the combination of op- 
tical and infrared data enables derivation of photometric red- 
shifts. The original plan for this program was to use the joint 
FLAMEX and NDWFS data sets to derive photometric red- 
shifts, and then detect galaxy overdensities within redshift 



slices. We developed a detection technique within this frame- 
work that is generically applicable to any photometric redshift 
catalog. The photometric redshift catalog is used to subdivide 
the data set into overlapping redshift bins, and we then employ 
a wavelet algorithm with bootstrap resampling to identify sta- 
tistically significant overdensities on a fixed physical scale (a 
complete description will be provided in Gonzalez et al., in 
prep.). 

Over the course of the survey this original plan has evolved 
in one significant regard. Joining forces with the IRAC Shal- 
low Survey (Eisenhardt et al. 2004), we are now using the 
combination of NDWFS, FLAMEX, and IRAC (3.6^m and 
4.5/im) photometry to derive more robust photometric red- 
shifts (Brodwin et al., in prep.) and probe more deeply than 
would be possible with either the NDWFS and FLAMEX or 
NDWFS and IRAC data sets alone. The IRAC data are deeper 
than FLAMEX, allowing us to push further down the cluster 
luminosity function, while the FLAMEX data significantly re- 
duce the scatter in photometric redshifts at z > 1 .4 (Brodwin 
et al., in prep.). We also now fold in the full redshift proba- 
bility distribution for each source, which minimizes the noise 
from galaxies with poorly constrained redshifts. Additionally, 
our joint team is using the same technique in a complemen- 
tary survey to detect clusters within the full 9.3 deg^ NDWFS 
Bootes field using only the NDWFS and IRAC data (Eisen- 
hardt et al., in prep.). 

In this paper we present one of the first spectroscopically 
confirmed z > 1 clusters detected with the combined NDWFS, 
FLAMEX, and IRAC data sets, CL1432.5H-3332.8 at z = 1 . 1 1 
(Figure im. To identify clusters within the FLAMEX region 
we perform a wavelet search using the photometric redshift 
probability distributions derived from the combined data sets. 



9 




Fig. 11. — The four image panels, each 3' X 3', show the core region of CL1432.5+3332.8 in Bw, L Ks, and IRAC 4.5 fim. The cluster is visible as an excess of 
red galaxies in the Kg and 4.5/xm bands, and in the A^,— band image we overlay circles to denote all galaxies with photometric redshifts 1 < z < 1.2. The centroid 
of the cluster detection is located at a = 14h32m29.2s, <5 = +33d32m48s (J2000), which is near the dense clump of red galaxies in the upper, central part of the 
images. 
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Fig. 12.— Spectrum of a galaxy at z = 1.110 in CL1432.5+3332.8. after 
smoothing with a 5 pixel boxcar filter. Several absorption features common 
to early-type galaxies are marked. 



with the wavelet scale size defined to detect structures with 
physical scales r « 200-800 kpc. Monte Carlo simulations 
are then used to quantify the significance of each detection. 
This search identifies the cluster CL1432. 5+3332. 8 as one of 
the three highest significance overdensities at z > 1, with an 
estimated redshift z = 1 .2 (FigurelTTl. 

Spectroscopic confirmation of this candidate was obtained 
on 3 June 2005 UT using LRIS on Keck 1 . We used a single 
slit mask for the observations which were carried out in clear 
conditions with 0.9 arcsec seeing. Five 1800s exposures were 
obtained, with the target galaxies shifted within the slitlets be- 
tween each exposure. T he data were re duced using standard 
long slit techniques (e.g. lStanford et al] ri997). The resulting 
redshifts include 8 galaxies at 1.10 < z < 1.12. Of these, five 
exhibit absorption line features indicative of old stellar pop- 
ulations (Figure [T2t. arguing that we are observing a cluster 
with an established early-type population. A histogram of the 
redshifts is shown in Figure [O] The resulting velocity dis- 
persion is (T = 920 ± 230 km s"', as calculated using the gap- 
per method (Beers, Flynn, & Gebhardt 1990), confirming that 
CL1432.5H-3332.8 is a massive galaxy cluster 

The cluster CL1432.5H-3332.8 is presented here as a proof 
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Fig. 13. — Redshift distribution from Keck observations of 
CL1432. 5+3332. 8. The solid histogram includes only galaxies with pho- 
tometric redshifts 0.95 < z < 1.25, while the open histogram includes all 
spectroscopic sources. The derived redshift and dispersion are < z >= 1.109 
and o- = 920 ± 230 kms"'. 



The survey area (7.1 deg^) is divided between the north- 
ern (Bootes) and southern (Cetus) NDWFS regions, with 
4.7 deg^ in Bootes and 2.4 deg^ in Cetus. We provide an ini- 
tial characterization of the data set based upon the 150,000 
sources detected (5a) within the Bootes field. The data in this 
region reach a median depth of Ks ^ 19.3 (6" aperture), with a 
positionally dependent limiting magnitude that we model us- 
ing extensive Monte Carlo simulations. Using the Bootes cat- 
alog in conjunction with optical NDWFS photometry we are 
able to perform a color-based star-galaxy separation that is 
robust for all sources except faint, blue galaxies at z 0. We 
demonstrate that the Bootes number counts are fully consis- 
tent with results from previous, smaller area surveys down to 
Ks = 19.6, but with vastly reduced statistical uncertainty. We 
discuss the utility of this data set for quantifying the number 
density and clustering properties of the ERO population. The 
FLAMEX Bootes catalog alone contains an order of magni- 
tude more EROs than any published studies. Finally, we de- 
scribe one of the first spectroscopically confirmed clusters at 
z > 1.1, which was identified using photometric redshifts de- 
rived from the joint FLAMEX, NDWFS, and IRAC Shallow 
Survey catalogs. 



of concept to demonstrate that we can successfully identify 
clusters at z > 1 using the combined Bootes data sets. A sin- 
gle cluster detection of course cannot constrain the fidelity or 
completeness of catalogs derived with this data set. These are 
issues that will be addressed in subs equent papers, and in par- 
ticular we refer the reader to Stan ford et aLl ( 120051) . Brodwin 
et al. (in prep.) and Eisenhardt et al. (in prep.) for fur- 
ther explorations of the photometric redshifts and cluster cat- 
alogs derived from the NDWFSh-FLAMEXh-IRAC and ND- 
WFSh-IRAC data sets. 

5.3. Public Catalogs 

The FLAMINGOS Extragalactic Survey received time 
through the NOAO Large Survey Program and pro- 
vides public access to the catalogs for the full FLAMEX 
region. The most recent version of the catalogs, ac- 
companied by a detailed description, is posted on the 
web (starting September 30, 2005) at the following mirrors: 
http : / / flamingos .astro.ufl.edu/ extragalact 
and ,http : / / www . noao . edu/ noao/noaodeep/] 

6. SUMMARY 

We present an overview of the FLAMINGOS Extragalac- 
tic Survey (FLAMEX), a near-infrared imaging survey con- 
ducted at the KPNO 2.1m with the FLAMINGOS instrument. 
FLAMEX is the first deep, wide-area NIR survey to image in 
both the / and K, filters, and is larger than any previous sur- 
veys to comparable depth. Located within the NOAO Deep- 
Wide Field Survey regions, FLAMEX is designed to com- 
plement the NDWFS optical photometry and other existing 
multiwavelength public data sets, with the intent that it be a 
resource to the astronomical community. We therefore make 
the full catalogs fro m thi s program publicly available via the 
web pages listed in ^5.31 
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